Scaling of dynamics in 2d semi-dilute polymer solutions 
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We study the viscoelasticity of surface polymer monolayers by measuring the dynamics of thermal 
concentration fluctuations with surface light scattering. For various systems of proteins and synthetic 
polymers we find a semi-dilute regime in which both the elastic and viscous components of the 
dilational modulus show power law dependencies on the concentration. Surprisingly there is a 
universal relationship between the exponents for the two components: the viscosity scales with 
a power double that of the elasticity. These results cannot be explained on the basis of theory 
developed for bulk systems, and a simple explanation for the singular 2d behavior is suggested. 
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Polymers are flexible long chain molecules, of outstand- 
ing importance in diverse fields from processing of ma- 
terials to biological activity. It is well known that it is 
possible to effectively confine some polymers to two di- 
mensions, for example by anchoring each monomer to the 
interface between immiscible fluids Q . Polymer dynam- 
ics in two dimensions (2d) has remained relatively unex- 
plored compared to bulk solutions, despite there being 
important examples both in life sciences and in technol- 
ogy where polymer molecules are confined to a plane. We 
consider the simplest situation where the chains lie flat 
and are free to move only in a 2d space. Such a surface 
layer exerts a lateral osmotic pressure EE = 70 — 7, which 
is the amount by which the surface tension 7 is reduced 
compared to the free interface tension 70. This pressure 
becomes significant above overlap of single chains, where 
it has a power-law dependence on the concentration that 
has been known for some time 0]. Like its counterpart 
in 3d bulk solutions, this equilibrium property is well 
explained within the picture of the polymer semi-dilute 
regime as an ensemble of independent "blobs" De- 
spite a few decades of experiments, very little is known 
instead about the physical nature of the divergence of vis- 
cosity and elasticity in 2d. These dynamical properties 
control important processes like foam drainage and 
stabilization and flow in emulsions when polymers are 
used as surfactants. The theoretical framework to de- 
scribe complex non-Newtonian flow and the dynamics of 
chains is well understood in 3d , but cannot be applied 
straightforwardly in 2d. How does the confinement to the 
surface, specifically a 2d effect, affect the dynamics of the 
system? 

This Letter presents a study of the dynamics of con- 
centration fluctuations in a wide range of polymer mono- 
layers, performed with surface dynamic light scattering. 
Polymers are spread onto an interface in dilute condi- 
tions, and cannot subsequently submerge and re-surface. 
Under these conditions polymer chains are unable to 
cross each other. This is the case for many monolayer 
systems studied in the literature. The technique used 
in this work, developed by Langevin and others, see the 



monograph [jj , is the only existing probe of thermal con- 
centration fluctuations in monolayers. Despite its use on 
polymer monolayers by various groups [8j, some issues 
concerning the data analysis have been resolved only re- 
cently, hence the data which is available for a quantitative 
comparison is very limited. 

Isolated polymer chains are usually modelled as ran- 
dom walks with potential interactions Q. The Flory 
exponent v relates the number of monomers N to the 
radius of gyration: R g ~ N v . When the concentra- 
tion in a polymer solution is increased so that individual 
chains are forced to overlap, the system enters a semi- 
dilute regime that lasts until the monomer fraction is 
very high. The semi-dilute region is important because 
for sufficiently long chains it covers a wide range of con- 
centrations. The equilibrium properties of polymers in 
this regime are given both in 2d and 3d by scaling laws, 
with exponents related to v. In particular, fluctuations 
of the density of monomers are correlated over a length £, 
defining a region known as a "blob" [3} . Inside the blob 
the chain has the same statistics as an isolated single 
chain, but different blobs are statistically independent. 
In 2d £ scales with the concentration V as: 

decreasing from a value of the order of R g at the over- 
lap concentration to the monomer size. Scaling of the 
osmotic pressure in a monolayer was first shown experi- 
mentally in ref. 0]: 

n ~ C 2 ~ r Ve *, where y eq = 2v/(2v - 1). (2) 

Figure ^ shows equilibrium data for some of the systems 
studied in this work. In 2d there are well known lim- 
iting regimes: The "good solvent", in which the chain 
behaves as a self-avoiding random walk, hence v = 3/4 
and y eq = 3, and the "6* conditions" where the effects of 
excluded volume balance the monomer-monomer prefer- 
ential attraction. Here the Flory exponent is predicted to 
be v = 4/7 and the power law exponent increases to 
y eq — 8, making this isotherm much steeper and compact. 
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FIG. 1: Surface pressure II as a function of area. Data 
are representative of three systems studied in this Letter: 
(o) PVAc, (V) /3-casein, (a) /3-lactoglobulin. Segments in- 
dicate the semi-dilute regime region where the equilibrium 
scaling exponents are determined. Inset: log plot of the equi- 
librium dilational modulus e eq against concentration. Lines 
are fitted power laws, showing the well known scaling of this 
equilibrium property in the semi-dilute regime. 



FIG. 2: (color online) The geometry of scattering from ther- 
mal surface roughness fluctuations, as probed by surface light 
scattering (SQELS). Typical values of the scattering vector q 
are of the order of 400cm -1 , corresponding to capillary wave 
frequencies lo of order 10 5 Hz. A is the frequency of the laser 
light. Details of the technique are given in ref. 



wave frequency to the wavelength has to be known: 



In contrast to the bulk, where the polymer chains can in- 
terpenetrate, in the planar 2d geometry the polymers are 
segregated and attain disk-like configurations |3| • Chains 
are corralled by their neighbors, as was observed with 
DNA molecules confined to a plane in ref. [Io| . 

In general, the response to a deformation in an 
isotropic 2d material is characterized by two elastic mod- 
uli: changes in area are controlled by the dilation modu- 
lus £ and changes in shape by the shear modulus G 
Polymer monolayers in the semi-dilute re gim e are fluid- 
like and the shear modulus is negligible 12]. In these 
conditions it is common practice to determine the equi- 
librium dilational modulus e eq (the 2d analog of the bulk 
modulus) from measurements of pressure as a function 
of = dH/dlnA. By definition, e eq has the same 

scaling properties as the osmotic pressure, described in 
Eq. |21 Power laws of function of concentration 

are shown in Fig. ^ The dilational modulus can also be 
measured dynamically, in an experiment where the sur- 
face area oscillates in time with a frequency to. Then 
the dynamic complex modulus e* — s' + itoe" is probed, 
where e' is the elastic component of the response modu- 
lus and e" is the dilational viscosity, e* can be accessed 
with the surface quasi-elastic light scattering (SQELS) 
technique, shown in Fig. [21 SQELS measures the time 
correlation function of light scattered from thermal sur- 
face roughness, which acts as a phase grating These 
out of plane fluctuations are underdamped waves with 
a frequency u>, and their motion is affected by the pres- 
ence of a surface film. To recover surface viscoelaticity, 
the surface wave dispersion relation D(to) that relates the 
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X] is the subphase Newtonian viscosity and p is the sub- 
phase density. Buzza 0] recently proved that the model 
Eq. |3| which is commonly found in the literature 0, El i 
is correct but that even under dynamical conditions the 
surface tension 7 should be considered as a real quantity, 
equal to the equilibrium static surface tension [23. As 
derived in ref. 15], the spectrum of light scattered by 
thermal roughness is given by: 
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SQELS data is fitted with Eq. [3| and Eq. H with only 
three physical parameters to be determined. Details of 
our experimental methods, including calibration, data 
analysis and limits of the technique, are described in 
ref. 0. 

In polymer systems, different models are appropriate 
depending on the lengthscales and timescales that are be- 
ing observed. The lengthscales q~ x probed with SQELS 
satisfy q£ < 1, so the semidilute solution is expected to 
be in a "macroscopic" regime where it behaves like a 
gel 0| • The gel's response is characterized by a rigidity 
modulus, e', scaling like the number of contacts, hence 
proportional to the osmotic pressure, and by a viscous 
dissipation modulus toe" describing the friction involving 
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FIG. 3: Log plot of the dilational elastic response modulus 
measured with SQELS as a function of the concentration, 
showing power-law scaling. Open symbols are the elastic (e') 
and filled symbols the viscous (we") components. Symbols 
are for different monolayers: (o) PVAc, (V) /3-casein, (a) (3- 
lactoglobulin, corresponding to data in Fig.0 Standard devi- 
ations obtained from fitting repeated SQELS correlation func- 
tions (not shown on each data point for clarity) are used to 
weigh the power-law interpolation (lines) . The exponents ob- 
tained from these fits are given in Table Eland are plotted in 
Fig- El The correlation between the different exponents, high- 
lighted in Fig. 01 represents the major result of this Letter. 



the monomers and the solvent. This is well established in 
three dimensional solutions, where both a fast and a slow 
relaxation mode can be probed simultaneously, for ex- 
ample with dynamic light scattering . The fast relax- 
ation can be described within the Rouse model and is re- 
lated to a cooperative diffusion timescale T coop 17]. This 
describes fluctuations that become faster as the concen- 
tration increases and the correlation length £ decreases. 
The slow relaxation is related to the self diffusion coef- 
ficient for reptation and is determined by the time T rept 
required for the chain to diffuse along its length. r rept 
is an increasing function of the concentration. The ap- 
proach of directly measuring the frequency spectrum of 
light scattered by concentration fluctuations is not pos- 
sible in monolayers because there is insufficient optical 
contrast, and dynamics can be probed only by indirect 
methods such as SQELS. 

Figure |31 shows the components of the monolayer vis- 
coclatic moduli, measured with SQELS for different poly- 
mer monolayers. Both the elastic and viscous compo- 
nents of e* exhibit a power law dependence on the con- 
centration: 

e' ~ r y =' and e" ~ T Vc " (5) 

A similar scaling behavior was recently reported by Mon- 
roy et al. 0] but has not been explained [24| . Results for 
all the systems considered in this work are summarized 
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Ref. 


Veq 


y e > 


Ve" 


PVAc 


T=45°C 


[16] 


2.0 


2.0 


5.9 


PVAc 


T=25°C 


[16] 


1.9 


2.4 


5.5 


PVAc 


T=6°C 


[16] 


1.7 


2.4 


5.5 


PVAc 


T=25°C 


[19] 


2.8 


2.8 


4.3 


/3-casein 


pH=5.3, 0.01M 


[201 


7.1 


5.6 


13.7 


/3-casein 


pH=7.2, 0.01M 


[20J 


6.0 


5.8 


11.9 


/3-casein 


pH=8.3, 0.001M 


[20] 


5.6 


6.0 


12.8 


/3-casein 


pH=8.3, 0.01M 


[20] 


5.3 


5.3 


10.8 


/3-casein 


pH=7.6, 1.1M 


[20] 


4.6 


4.3 


9.3 


/3-lg 


pH=6.0, 0.02M 


[22] 


8.3 


5.4 


15.1 


/3-lg 


pH=8.3, 0.02M 


[22] 


5.3 


6.2 


13.3 


/9-lg 


pH=5.9, 0.1M 


[22] 


8.8 


7.6 


13.2 


P4HS 


P H=2, T=25°C 


fl9l 


8.1 


9.0 


14.3 



TABLE I: Summary of static and dynamic scaling exponents 
for monolayers studied in the literature. Values of y eq are ob- 
tained from measurements of the equilibrium surface pressure 
with Langmuir trough methods, while values of y e i and y E n 
are obtained from SQELS data. The SQELS data of Cicuta 
and Hopkinson has been fitted following the method described 
in ref. The references specified give experimental details 
for each monolayer. 

in Table d Monolayers of very different compositions, 
comprising both synthetic polymers and proteins, have 
been studied, to cover the widest possible range of values 
of v. The temperature, ionic strength and pH of the liq- 
uid subphase are all controlled, as they affect the chain 
configuration at the surface. 

This data contains a wealth of information. Looking at 
Fig. it can be clearly seen that for each monolayer the 
viscosity scales with a higher power than the elastic mod- 
ulus. On general dimensionality terms, the viscosity can 
be regarded as the product of the modulus driving the 
relaxation and a characteristic time U s" ~ e' ■ r. From 
this it immediately follows that the timescale r of these 
relaxations becomes longer the higher the concentration, 
hence a slow mode is being probed. The character of this 
mode can be read from Figure 01 where the power law 
exponents for e* defined in Eq. are plotted against the 
equilibrium exponent. It is not surprising to find that 
Ve 1 — Veq, since the thermal fluctuations probed with 
SQELS fall within linear response theory. What is of in- 
terest is the apparently universally valid relation between 
the exponents for the elastic and the viscous components 
of e*: y £ » — 2y £ /. It finally follows that the timescale 
of these fluctuations, which we shall call T2d, scales like 

r 2d ~ r 2 - 

We are only able to speculate on the physical origin 
of this timescale. A simple explanation is that the slow- 
ing down of T2d as the concentration increases is due to 
friction between the increased number of statistically in- 
dependent blobs. The number of contacts in between 
blobs is proportional to £~ 2 , giving: 

e" „ r 4 ~ T 2 ^' (6) 
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FIG. 4: The scaling exponents describing the power law 
dependence on the concentration of the dilational elastic- 
ity (open symbols) and viscosity (filled symbols) are plot- 
ted against the exponent for the equilibrium bulk modulus. 
(•): PVAc, (T ): 3 -casein, (A): /3-lactoglobulin, (■): data of 
Monroy et al. |1S1| - This figure includes results from the data 
shown in Fig.[3]as well as the monolayers described in Table[l] 
The lines have slopes of 1 and 2. It is clear from the data that 
the exponent for the viscosity is approximately twice the ex- 
ponent of the elastic modulus. This result is different from 
analogous measurements on bulk polymer solutions, indicat- 
ing a specific 2d dynamics. 

which describes the data of Fig.0]very well. An analogy 
can be drawn between the diverging viscosity of the close- 
packed arrangement of blobs in the semi-dilute regime 
and the case of diverging viscosity of spheres at high 
packing density \2l\ , both being determined by the num- 
ber of contacts. 

A question that arises is why this mode is not seen in 
3d. In bulk solutions in a good solvent a dynamical mode 
analogous to that just described for 2d would have a char- 
acteristic relaxation time scaling with the concentration 
(j) like t ~ £ -3 ~ 9 / 4 . This is a higher power of the 
concentration compared to reptation (r revt ~ 1 ' 5 )0|j m 
agreement with the well known result that in 3d solu- 
tions the fluctuations relax via the self diffusion of the 
polymer chain and not the mechanism outlined above. 
In 2d relaxation by reptation is hindered by an effective 
confinement provided by neighboring chains. 

In summary, it has been shown that the thermal con- 
centration fluctuations in the semi-dilute regime of poly- 
mer monolayers reported in the present paper and in 
ref. 01 have a slow mode of decay with a timescale that 
had not been previously considered and that describes a 
specifically 2d process. 

We thank E.M.Terentjev for very useful comments and 
discussions. 
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